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CjH' ABSTRACT 

We present a library of 6410 synthetic spectra with resolution A/AA = 250 000 based 
on the revised Kurucz 1993 model atmospheres. The library covers the wavelength 
range 3000 — 10 000 A with 54 values of effective temperature in the range 5250 — 50 000 
K, 11 values of log surface gravity between 0.0 and 5.0 and 19 metallicities in the range 
—5.0 to 1.0. We find that, with a few caveats, the library compares well with both the 
original 20 A Kurucz spectra and also with observed spectra. The library is intended 
for use in population synthesis and physical parameterisation of stellar spectra. We 
assess the suitability of the library for these tasks. 
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Progress in population synthesis and automatic classifica- 
tion of stellar spectra has been limited by the spectral res- 
olution of the available synthetic stellar spectra. The ex- 
isting synthetic libraries are not at high enough resolu- 
tion to be useful f or classifying stars fro m recent surveys 
such as the SDSS JStoughton et alJl2002l) CA/AA ~ 1 800), 
or futu re surveys such as RAVE |Steinmet3 [2002) and 
GAIA iLindeeren fc Perrvmanl Il996l) . Most classification 
techniques smooth the observed spectra to the resolution 
of the synthetic spectra. This means much of the detailed 
information in the observed spectra is lost, which may re- 
duce the quality of the cl assifications. Population synthesis 
packages s uch as pegase flFiocfc Rocca- Volmerangd Il997ft 
or gissel feruzual A. fc Charlodll993l) use a grid of stel- 
lar spectra to generate galaxy spectra. The resulting galaxy 
spectra are limited by the resolution of the input stellar 
spectra. To study the high resolution features, it is neces- 
sary to have a grid of observed or synthetic stellar spectra 
which match the resolution of the observed galaxy spectra. 
For example, galaxy spectra synthesised from 20 A spectra 
cannot be used to measure standard line indices like the Lick 
indices. 

Perhaps the most widely used library of synthetic spec- 
tra are the flu x distributions from the Kurucz ATLAS9 model 
atmospheres dKurucz) ll993aD . It is important to note that 
while these are usually referred to as spectra, they are 
flux distributions predicted directly from the model atmo- 
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spheres, rather than spectra generated by a spectral syn- 
thesis program. The Kurucz atmospheres have several dis- 
advantages w hich are discussed in various sources such as 
iKuruczl il992T) . However, one of their advantages is the wide 
range of parameters they cover, which is important for gen- 
erating a grid of stellar spectra for population synthesis. 
The need for higher resolution spectra has been recognised 
for some time but, because of the immense computational 
expense involved, the synthesis of the spectra has been lim- 
ited to partial wavelength ranges and specific regions of pa- 
rameter space. Several groups have generated libraries of 
spectra from the Kurucz mode l atmo spheres. For example, 
IChavez. Malagnini. fc Morossil (119971) provide a set of 711 
Kurucz spectra at A/AA = 250 000 in the w avelength region 
4850 - 5400 A and lCastelli fc Munaril (l200lt) have generated 
a set of 698 Kurucz spectra at A/AA = 20 000, in the wave- 
length region 7650 - 8750 A for use with GAIA spectra. 
iGonzalez Delgado fc Leithererl Jl99Sl) have generated syn- 
thetic spectra in very small spectra regions necessary for a 
particular application. They created a grid of synthetic pro- 
files of stellar H Balmer and Hel lines at AA = 0.3 A for the 
purposes of evolutionary synthesis. 

There are also several libraries of observed spectra now 
available at much higher resolution. Fo r example, the elodie 
database JPrugniel fc SoubiranllioOlh . consists of 709 stars 
observed in the wavelength range 41 00 — 6800 A with a res- 
olution of A/AA ~ 42 000. STELIB ihe Borgne et alJl2003h 
provides spectra for 249 stars observed in the wavelength 
range 3200 - 9500 A with a resolution of A/AA ~ 2000. The 
observed libraries are crucial for evaluating the accuracy of 
synthetic spectra and can also be used directly for popu- 
lat ion synthesis and classifi cation. STELIB has been used 
bv lKauffmann et al.l <l2003l) with a new version of GISSEL, 
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and elodie has been used to assign physical parameters to 
stars observed by the SDSS. However, a limitation of the 
observed spectral libraries is that they do not cover the full 
range in parameter space needed for galaxy population syn- 
thesis. Complete coverage of the parameter space is even 
more important for stellar spectral classification. The most 
successful approaches t o classification hav e used methods 
from machine learning jBailer-Jonesll200lL In these meth- 
ods, the distribution of spectra in the training set has a 
direct impact on the accuracy of the classification assigned 
to new spectra. 

There are several efforts to generate hi gher resolution 
Kuruc z spectra that are currently in progress. iBertone et alJ 
(2002) have generated a grid of 832 spectra at A/AA = 
500 000 over the wavelength range 3500 - 7000 A. They in- 
tend to extend the wavel ength range down to 850 A at areso- 
lution of A/AA = 50 OOO. IZwitter. Munari. fc Castellil (120021) 
are in the process of generating a grid of Kurucz spectra at 
A/AA = 20 000 over the wavelength range 2500 - 10 500 A 
for use in radial velocity correction work. 

We have generated a larger library of 6410 spectra from 
the Kurucz model atmospheres. Previously these spectra 
were only available either at much lower resolution (20 A) or 
over small wavelength ranges. Our spectra were generated 
from ATLAS9 model atmospheres, using John Lester's Unix 
version of the SYNTHE spectral synthesis package (Lester 
2002, private communication). We have modified this pack- 
age to improve the efficiency of the code, making it possible 
to generate the complete range of Kurucz spectra in a rea- 
sonable time. 

This paper compares our higher resolution spectra with 
the original 20 A Kurucz spectra and the STELIB library of 
observed spectra. In Section[5]we describe the main charac- 
teristics of the new library of spectra. In Section El we com- 
pare th e spectra with the 20 A Kurucz spectra from lKuruczl 
lll993al) . Finally, in Section [I] we compare the spectra with 
observed spectra from the STELIB library. We will make 
this library of spectra available for general use on request. 



2 GENERATING THE KURUCZ SPECTRA 

The library presented here has been created using the 
updated versio ns of the ATLAS9 model atmospheres from 
IKuruczl ljl993al) . These are available from Kurucz 's website 
(labelled .dat). Kurucz advises that the iKuruczl (ll993aT) 
models (labelled .datCD) should not be used as they have a 
discontinuity in the fluxes and colours as a function of T e a 
and log(<?) that was corrected for the revised version. 

The models assume plane parallel homogeneous lay- 
ers in steady-state, local thermal equilibrium. A microtur- 
bulance velocity of 2 k ms -1 and a mixing-length value o f 
£/H p = 1.25 are used. ICastelli. Gratton. fc Kuruczl (Il997l) 
give a detailed discussion of whether the mixing length the- 
ory for convection is dealt with adequately in the stan- 
dard Kurucz atmospheres. They have calculated an alterna- 
tive set of atmospheres in which the Kurucz "overshooting" 
approximation is not used (NOVER models). These mod- 
els have been shown to predict more accurat e observable 
properties (eg. colours) for some atmospheres (j|H eiter et alJ 
Eflfla ISmallev fc Kur,"kHl99fl ISma.llev et a.l]l2002hL How- 
ever, since they are currently available only for some metal- 



licity values (-2.5, -2.0, -1.5, -1.0, -0.5, 0.0, +0.5 dex) we 
decided not to use these atmospheres for population synthe- 
sis. We have generated a small subset of the NOVER models 
for comparison with the STELIB spectra. 

The adopted atomic line lists for all our sp ectra are 
LOWLINES and NLTELINES from lKuruczl l)l994t) and the 
adopted molecular line list is DIATOMIC from IKuruczl 
dl993bl) . We have included both predicted and measured 
lines. The inclusion of the predicted lines is necessary to re- 
produce accurate flux distributions. However, it does mean 
that care should be taken when using the high resolution 
spectra, as the properties of individual lines may not be as 
accurate as if Kurucz's more up to date linelists (which do 
not include predicted lines) were used. The molecular lines 
have been included for spectra with T a a ^ 7000 K. 

The original programs making up the Kurucz package 
were written in F0RTRAN77 for the VAX. John Lester has writ- 
ten a Unix version based on the original Kurucz code, in 
F0RTRAN77 and a more recent Unix version in Fortran90 
(Lester 2002, private communication). However, this code 
could only be used to generate small sections of a spec- 
trum at a time because it required massive quantities of disk 
space. The complexities in generating spectra with SYNTHE 
has restricted researchers to creating either a small number 
of spectra, or a large number over a very small wavelength 
range. We have modified the Lester code to dramatically 
reduce the disk usage of SYNTHE which makes it feasible 
to generate large spectral ranges on standard hardware in 
a relatively short period of time. We have not yet gener- 
ated spectra for values of T E ff ^ 5000 K. These spectra 
require TiO lines to be included which makes the program 
significantly slower, and requires more disk space and we are 
investigating whether further optimisations may be made. 

One thing to note is the subtle difference in terminology 
between 'spectra' and 'flux distributions' which are some- 
times used interchangeably. The Kurucz 20 A flux distribu- 
tions were predicted directly from the Kurucz model atmo- 
spheres. The ATLAS9 model atmosphere program sums up 
the opacity in broad 10 or 20 A bins creating a low res- 
olution spectrum. The Kurucz spectra that we have gen- 
erated have been calculated from the model atmospheres, 
using the SYNTHE software. These spectra are generated at 
extremely high resolution (ideally as high as the compu- 
tational resources allow) and then rebinned to the required 
resolution. The newest versions of ATLAS9 will combine both 
the model atmosphere generation and the spectral synthesis 
into one process. 

Of the 7216 model atmospheres available in Kurucz's 
standard distribution, we have generated a grid of 6410 
spectra, which excludes the lowest temperature spectra. The 
spectra cover the wavelength range 3000 — 10 000 A which 
was chosen to be useful for comparisons with SDSS spec- 
tra. The spectra were generated at a resolution of A/AA = 
250 000. 



3 COMPARISONS WITH THE 20A KURUCZ 
SPECTRA 

In this section we compare our library of spectra with the 
revised version of the 20 A Kurucz flux distributions. There 
have been several changes to the software and data between 
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Figure 1. The distribution of our Kurucz spectra in log(T c fj ) - 
log(g) space (x). Also plotted is the distribution of the STELIB 
spectra with T eB > 5000 K (+). 



the original release of the 20 A Kurucz flux distributions, 
and the present. Also, the flux distributions are predicted 
directly from the model atmosphere code, so it is not pos- 
sible to compare them directly with the spectra generated 
by spectral synthesis. Because of this, the spectra we have 
generated are not expected to match the flux distributions 
exactly, but it is an important check to make sure there was 
broad agreement between the spectra. 



3.1 Direct comparisons 

We compare the higher resolution Kurucz spectra with the 
20 A flux distributions. Figure |5] shows several examples of 
these comparisons. For the purposes of this comparison, we 
re-binned the new spectra to 20 A using a simple top-hat 
function. There are several differences between the two sets 
of spectra: 

(i) In all the spectra there is an extra dip in the new spec- 
tra at around 3700 A. This is due to inadequate treatment 
of the Balmer jump. In the model atmosphere code, the way 
the Balmer jump is dealt with has been modified in order 
to remove the dip (Kurucz, private communication). 

(ii) In some of the lower temperature spectra, the shape 
of the Ca H/K doublet is not identical. This is largely due to 
the fact that for the model atmospheres, the Stark broaden- 
ing was artificially increased to make up for missing lines in 
computing the distribution function (Kurucz, private com- 
munication). This was changed back to generate the syn- 
thetic spectra, and also for more recent flux distributions. 

(iii) In lower temperature spectra (T e g ^ 7000 K) there is 
a systematic difference in flux in the G-band. This is likely 
to be due to changes in the molecular line lists since the flux 
distributions were produced. 

(iv) Most of the strong lines have slightly different depths. 
These small peaks are likely to be unresolved line cores, 
caused by the fact that the 20 A flux distributions are under- 
sampled, whereas our new synthetic spectra are generated 
at high resolution. 



Probably the most significant of these points for our ap- 
plications is the difference in the G-band flux. This affects 
individual lines, and the U band magnitude of the spectra. 
In the following sections we quantify this. 

3.2 Colours 

Comparing the UBVRI magnitudes is important for check- 
ing that the broadband properties of the spectra are repro- 
duced. We have used th e Johnson-Cousins bandpasses, as 
defined bv lBesselll Jl990|) . Fi gure |21 shows the offset in mag- 
nitudes calculated for the new spectra and the 20 A flux dis- 
tributions. The RMS errors for the colours for the two sets 
of spectra are: AU = 0.05 mag, AB = 0.01 mag, AV = 0.01 
mag, AT? = 0.01 mag and AI = 0.01 mag. The B, V, R and 
I magnitudes are in good agreement. However, the differ- 
ences are greater for the U magnitudes at lower tempera- 
tures (T e a ^ 7000 K). For higher temperature spectra the 
U magnitudes are also in good agreement. The U band dis- 
crepancy is probably a consequence of the differences in the 
spectra noted in Section f3. II 

The question is how important these magnitude offsets 
are for the applications we are interested in. A difference in 
measured colour is lar gely degenerate with a differen ce in 
T e ff. For example. iLeieune. Cuisinier. fc Buserl (1997) find 
that the colour vs. temperature relations for Kurucz spectra 
do not match those derived empirically. To fix the prob- 
lem, they developed a method for adjusting the shape of 
the continuum (in other words change the magnitudes) for 
spectra of a certain T ft . However, an alternative approach 
is to leave the shape of the spectrum unchanged and re- 
assign a T e ff to the spectrum. For example, a solar metal- 
licity spectrum with T c s = 7500 and log(g) — 1.0 has a 
U magnitude offset of 0.05. This difference in colour be- 
tween our new spectrum and the Kurucz flux distribution 
corresponds to a change of T e s by less than 30 K. Since 
quoted errors for temperatures on observ ed spectra a r e typ - 
ic ally around 40 — 100 K (see for examplelKatz et alJ ill 998f) 
or lAlonso. Arribas. &i Martinez- Roger! (ll99STl ^l . the spectra 
agree to within usual measurement errors. 



4 COMPARISONS WITH OBSERVED 
SPECTRA 

Having shown that the new spectra agree reasonably well 
with the 20 A flux distributions, the next step is to evaluate 
them against observed spectra. This serves two purposes. 
Firstly, there are some differences between the two sets of 
Kurucz spectra, especially in the blue. However, if these dif- 
ferences are smaller than the differences between one of the 
Kurucz spectra and an observed spectrum with the same 
physical parameters, then they may be neglected. Secondly, 
we want to evaluate what kind of applications the synthetic 
spectra can be used for. Our aim is to demonstrate that our 
new Kurucz spectra are as suitable for use in population syn- 
thesis as the commonly used 20 A Kurucz flux distributions, 
but with the advantage of higher resolution. There are al- 
ready various comparisons of Kurucz spectra with observed 
spectra in the literature. For example, the c olour— T e g cor- 
relatio n for the spectra has been analysed bv lLeieune et all 
Jl997f) and the 20 A spectra have been compared with 
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Figure 2. Comparison between a sample of the 20 A Kurucz flux distributions (bold line) and the newly generated Kurucz spectra 
(thin line) (box-car smoothed to AA = 20A). These spectra cover a range of physical parameters: the title of each plot gives the values 
for T e ff, log(g) and [Fe/H]. The x-axis is the wavelength in A and the j/-axis is flux, F(X), in arbitrary units. Underneath each pair of 
spectra is AF = F(A) new - F(A) old . 
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Figure 3. Magnitude offsets between our Kurucz spectra and the 20 A flux distributions. The x-axis is logio(T e f f) and the y-axis shows 
the difference (in magnitudes) between the Kurucz spectra we have generated and the 20 A flux distributions. 



obser ved spectra bvlStra izys, Liubertas. fc Valia ugal jl997l) 
and IStraizvs. Lazauskaite. fc Valiaueal <|2002T) . We here 
make our own comparisons to focus on properties such as 
the Lick indices which are relevant to the specific applica- 
tions we are interested in. 

We have compared the sp ectra to those in the STELIB 
library l)Le Borgne et, all2 003l. which consists of 250 spectra 
in the visible range (3200 — 9500 A) with a spectral resolution 
of A A < 3 A and covering a wide range in parameter space. 
As well as a direct comparison of the SEDs, we compare 
the measured colours and Lick indices to confirm that the 
spectra meet out requirements. 

Many of the STELIB stars do not have values for all 
three physical parameters (T e g, log(<?) and [Fe/H]), have 
large section of the spectrum missing, or do not lie close to 
a point on our grid of spectra. Hence we have selected a sub- 
set of 125 of the spe ctra that have tabulated values for T e g, 
log(g) and [Fe/H] in lLe Borgne et al.l(|2003h . T cff > 5000/^, 
and which have a good x 2 match to the corresponding Ku- 
rucz spectrum. The values for m ost of the physical pa- 
ramet ers we re obtained from the _JCa^rel^e^trobelet^1 
<1997l) and ICavrel de Strobel. Soubiran. fc Ralitd (l200lD 
catalogue s. However, some were obta ined from the elodie 
database iPrueniel fc SoubiranllioOlh and several were cal- 
culated using the tgmet method llKatz et all Il998fl . We 
corrected the spe ctra for interstellar red dening using the 
values of A v from iLe Borgne et all (l200St) and the idl AS- 
trolib routine ccm_dered.pro. The STELIB spectra have 
not been corrected for atmospheric absorption, so they have 
strong absorption features (most noticeably in the Oxygen 
A-band around 7600A) which are not present in the syn- 
thetic spectra. 



4.1 Colours 



We have calculated colours for all of the STELIB spectra and 
their closest Kurucz matches. Since the wavelength range of 
the STELIB spectra is 3200 - 9500 A, we have calculated 
a modified U band magnitude (denoted U*), which is trun- 
cated at 3200 A for both the STELIB and Kurucz spectra. 
Figure [I] shows the magnitude offsets between the STELIB 
spectra and their Kurucz matches. The RMS error between 
the two sets of spectra are: AU* = 0.14 mag, AB = 0.05 
mag, AV = 0.02 mag, AR = 0.03 mag and AI = 0.04 mag. 
For comparison, we also calculated the magnitude offsets for 
the STELIB spectra and their closest Kurucz matches, using 
the 20 A flux distributions. The RMS error between these 
two sets of spectra are: AU* = 0.11 mag, AB = 0.05 mag, 
AV = 0.02 mag, AT? = 0.03 mag and AI = 0.04 mag. This 
demonstrates that, with the exception of U* , the scatter in 
the magnitudes is no worse for our new spectra than for the 
20 A spectra that are in standard use. 

Some differences are expected due to the mismatch be- 
tween the STELIB temperatures and those of the closest Ku- 
rucz match - many of the spectra have a difference in tem- 
perature of about 100 K, corresponding to a difference of a 
few hundredths in magnitude. However, the important point 
for our work is that our new spectra are as well matched 
to the observed spectra as the 20 A spectra that are often 
used for population synthesis. Figure shows various colour 
comparisons. There is some scatter in each of the plots and 
a slight systematic offset in R — I. Again however, this is 
no worse than when comparing the 20 A spectra with the 
STELIB spectra. 



6 T. Murphy and A. Meiksin 



0.4 
0.2 
0.0 
-0.2 
-0.4 



4 
J* 



0.4 r 
0.2 7 



-0.4 . 



+ + 



, 0.0: ^#*Mf*^4M^ ♦ + V 

-0.2 r + 



2 4 6 8 10 

B 



> 
i 




0.4 r 
0.2 r 



0.0 r +-H-+44 + #-^4f):+4%%# + +f+ + +H 



-0.2 r 
-0.4 r 



4 6 8 

R 



Figure 4. Magnitude offsets between the STELIB spectra and the closest Kurucz spectrum match from the grid. 



4.2 Line indices 



ing random errors are given by: 



We have calculated Lick indices for each of the 125 STELIB 
spectra in our sample, and for the matching Kurucz spectra. 
There are two types of index, atomic (/„) and molecular (I m ) 
which we have calculated from the standard formulae: 
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where A c i and A C 2 are the limits of the central bandpass 
denning the index (in A), S(X) is the object spectrum and 
C(A) is the linearly interpolated pseudo-continuum, denned 
by 
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and Xbi, Xb2, X r i and X r 2 are the limits of the blue and red 
continuum bands. The bandpass definitions for the Lick in- 
di ces and the red and blu e pseudo-continua are those defined 
in lWorthev et, alJ l)l994|) . These are given in Tabled 

A method for random err or estimation i n line -strength 
indices is outlined in detail bv lCardiel et al.l l|l99£T) . The re- 
sulting equations are based on a full analysis of the error 
propagation throughout the calculation process. The result- 
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for the molecular indices, where 
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Note that lCardiel et all lll998l) assume the size of each pixel 
dX (0 in their notation) is fixed, and so take it out of all 
the summations. However, in general (and in our case), dX 
is not fixed, and so we keep it inside the summation signs. 

Figure |5] shows the correlation between the Lick indices 
calculated from the STELIB spectra and the Kurucz spec- 
tra. We have not calculated the two TiO indices since TiO 
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Figure 5. Top Left: U* - B vs. B - V for STELIB spectra (+) and their Kurucz matches (o). Top Right: U* - B colours for STELIB 
spectra (x-axis) and their Kurucz matches (j/-axis). The solid line is the line "x = y" for comparison. Bottom Left: Comparison of B — V 
colours. Bottom Right: Comparison of R — I colours. 
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125 


4083 
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- 5285 


650 
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150 
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150 
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375 
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000 
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875 
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625 
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625 


- 5796 


625 
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375 
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375 
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875 
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375 
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625 
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625 


5922 


125 
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125 
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Table 1. Bandpass definations for the Lick indices calculated here. Each index has a type A or M signifying atomic and molecular indices, 
respectively. Atomic indices are measured in A and molecular indices are measured in magnitudes. Note that these are air wavelengths. 
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Figure 6. Each plot shows the correlation between the Lick index calculated on the STELIB spectrum (ic-axis) and the closest matching 
Kurucz spectrum from our grid (j/-axis). The solid line is the line u x = y" for comparison. We have not included the TiO indices for 
reasons discussed in the text. The RMS error for each index is given in Table HI 
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Name Mean value Av. Abs. Err. Med. Abs. Err. Type 

M 
M 
A 
A 
A 
A 
A 
A 
A 
A 
M 
M 
A 
A 
A 
A 
A 
A 
M 



Table 2. Statistical comparison between the STELIB Lick in- 
dices and those of their Kurucz matches. Column 2 shows the 
mean value of the STELIB index. Column 3 shows the average 
absolute offset between the Kurucz and STELIB indices. Col- 
umn 4 shows the median absolute offset between the Kurucz and 
STELIB indices. 



lines are switched off in these spectra. The average absolute 
offsets (ie: Yl /fc ~ Js ) for the correlations are given in Ta- 
ble |^] It should be noted that the Kurucz spectrum paired 
with a given STELIB spectrum is not necessarily the best 
possible spectrum that could be synthesised to match the 
STELIB spectrum, but simply the spectrum from our grid 
with the closest match in physical parameters. This means 
it does not take into account any particular properties of 
the STELIB spectrum that could in principle be modelled 
and which may affect the line indices. However, the com- 
parison does give an indication of the level of mismatch to 
be expected when comparing observed stellar spectra with a 
grid of synthetic spectra like ours. Such mismatches are an 
inevitable consequence of any automated comparisons be- 
tween a grid of model spectra and the vast number of spec- 
tra being made available from large surveys like the SDSS. 

We have calculated Lick indices and their associated er- 
rors f or a set of 9473 SDSS galaxies from the Early Data Re- 
lease JStoughton et alJl2002fl . This allows us to compare the 
scatter between the Lick index values for the STELIB and 
Kurucz spectra, with the expected accuracy of the indices 
measured from observed galaxy spectra. We have found that 
the random error (as calculated by Equations [§] and asso- 
ciated with each index measured in the SDSS spectra, is of 
the same order of magnitude as the average offset between 
the Kurucz and STELIB index measurement. This suggests 
that the indices measured from galaxy spectra synthesized 
from the Kurucz spectra will have at least comparable accu- 
racy to those measured with SDSS, and so our grid of stellar 
spectra is well matched to the SDSS galaxies. 



5 DISCUSSION 

We have generated a grid of theoretical spectra from the 
Kurucz model atmospheres. Since the intended use of these 
spectra is in population synthesis and stellar classification, 
we have made several comparisons to check the validity of 
using the spectra for these purposes. The broadband prop- 
erties of the spectra compare well with observed spectra, as 
do the line index measurements. The comparisons do not 
guarantee the accuracy of the Kurucz spectra as models of 
observed stellar spectra. Rather, they demonstrate that our 
new, high resolution Kurucz spectra are as good for use in 
population synthesis as the commonly used library of 20 A 
spectra. The advantage of these spectra over the previously 
available Kurucz spectra is that they allow the modelling of 
spectral line features such as the Lick indices. 

When using these spectra at high resolutions, it should 
be recognised that the spectra were generated using line lists 
that include 'predicted' lines. This is necessary to reproduce 
the broadband colours of the spectra accurately. However, it 
does mean that many of the individual lines present at high 
resolutions do not have measured properties. Also, the line 
lists that we used (such as LOWLINES) are known to have 
problems with the values for specific lines. As this mostly 
involves weak lines, it should not present any difficulties for 
population synthesis at the resolution of the SDSS. More 
accurate properties for these lines could be obtained using 
alternative values, f or example from th e Vienna Atomic Line 
Database (VALD) llKupka et al.ll2000l) . 

Also, the 20 A Ku rucz flux distribution s are often used 
with t he corrections of ILeieune et alJ (I1997T) ; ILeieune et at.l 
( 1998) applied. The corrections are an attempt to calibrate 
the spectra by comparing the synthetic model colours with 
empirical stellar colours. We investigated applying the same 
corrections to our high resolution spectra, but found that 
artifacts were introduced when using this technique directly. 
The Lejeune corrections mainly affect the lower temperature 
spectra, which we have not included in this library. For most 
of our spectra, the corrections are negligibly small. 

We intend to extend the library to the lower tempera- 
ture models in which TiO lines become important and also 
generate spectra for the NOVER models as these are more 
accurate. We are in the process of building a large library 
for popu lation synthesi s - a h igher resolution version of that 
done by ILeieune et al l (^)97|) - which will also incorporat e 
the NextGen models lllauschildt. Allard. fc Baronlll999ah . 
jHauschildt et al.lll999bf) . 
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